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Abstract 

The adiabatic motion of a charged, spinning, quantum particle in a two - di- 
mensional magnetic field is studied. A suitable set of operators generalizing the 
cinematical momenta and the guiding center operators of a particle moving in a ho- 
mogeneous magnetic field is constructed. This allows us to separate the two degrees 
of freedom of the system into a fast and a slow one, in the classical limit, the rapid 
rotation of the particle around the guiding center and the slow guiding center drift. 
In terms of these operators the Hamiltonian of the system rewrites as a power series 

in the magnetic length Ib = \f§§ and the fast and slow dynamics separates. The 
effective guiding center Hamiltonian is obtained to the second order in the adiabatic 
parameter Ib and reproduces correctly the classical limit. 



short title: Adiabatic motion in a 2-dimensional magnetic field 
PACS: 03.65.-w, 02.90.+p 



1 Introduction 



The study of the adiabatic motion of charged particles in slowly-varying magnetic fields 
is a very well developed subject of classical mechanics Q. It is of utmost importance 
in plasma physics, in the study of particle motion in the geomagnetic field as well as in 
astrophysics. For these reasons several efforts have been devoted over the years in order 
to gain a deeper understanding of the formal structure underling the problem. In this 
manner the original method of directly subjecting the equations of motion to a standard 
perturbative treatment |]J, has been supported by the works of C. S. Gardner 0, H. E. 
Mynich || and R. G. Littlejohn [||, |5|, |6| developing Hamiltonian theory and yielding a 
systematic procedure allowing the construction of adiabatic invariants. While Gardner 
and Mynich based their approach on field lines coordinates and canonical transformations, 
Littlejohn proceeded by investigating the Hamiltonian structure of the system and by 
phase space (symplectic) geometrical techniques. In this paper we deal with the problem of 
the adiabatic motion of charged, spinning, quantum particle in a slowly- varying magnetic 
field. To this task we develop an algebraic technique which very much reminds the one 
used by Littlejohn in classical mechanics. Nevertheless we make use of a different set of 
variables and do not make explicit use of symplectic geometry. We restrict our attention to 
particles moving in a plane under influence of a time-independent perpendicular magnetic 
field. This is not only an expedient in order to illustrate our technique in a simple case. 
Although the planarity constraint may seems to be artificial in the classical context it is 
not in the quantum-mechanical one. Devices producing the confinement of electrons on 
a plane are in fact widely studied in physics, as example in the quantum Hall effect 0. 
When these systems interacts with arbitrary inhomogeneous magnetic fields a realistic 
analysis indicates the effective dynamics to be described by a Hamiltonian proportional 
to the two-dimensional Laplace operator in minimal coupling with an effective magnetic 
field represented by the normal component of the original magnetic field M. This is 
the kind of system with which we are going to deal with. Our perturbative analysis is 
shown to be consistent with some general results concerning the ground state of quantum 
particles moving in an arbitrary two-dimensional magnetic fields |]9], 10 . 



Before facing the quantum problem it is better to briefly recall the main features of the 
classical one. We start therefore by considering a classical particle of mass m and charge 
e moving in a plane under the influence of a normal uniform magnetic field of magnitude 
B. Let x = (x 1 ,^ 2 ) denote the position of the particle and v = (v 1 ^ 2 ) its velocity. As 
well known the particle follows a circular orbit of radius r B = ^§\v\ which center remains 
motionless. The frequency of the motion is the cyclotron frequency ujb = —• The center 
of the orbit is called guiding center and its coordinates X = (X l ,X 2 ) are related to 
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Figure 1: Adiabatic motion of a charged particle in a two-dimensional magnetic field. 



the position and the velocity of the particle by the relation X % = x % + ^§£ lk Vktf\ The 
stronger the magnetic field the smaller is the radius of the orbit and faster the particle. 
By increasing the magnitude B of the magnetic field, hence, the classical particle explores 
a portion of the plane which becomes smaller and smaller. This gives a reason of the why 
in classical mechanics the motion in a very strong magnetic field may be studied along the 
same lines of that in weakly inhomogeneous one. This fact is not immediate in quantum 
mechanics since even in a uniform magnetic field the wavefunctions are not localized by 
increasing the magnetic field. 

Let us consider now the case of a weakly inhomogeneous magnetic field of magnitude 
B(x). What happens is that the particle keeps on rapidly rotating around the guiding 
center X but, this time, the guiding center no longer remains motionless. It begins to 
slowly drift in the plane. The cyclotron frequency becomes, in addition, a function of the 
position lob(x) ~ eB ( x ) _ yy e can therefore distinguish two way of move of the system: 
the very fast rotation around the guiding center and the slow drift of the guiding center 
in the plane. From the point of view of the Hamiltonian theory it is not immediate to 
understand how to deal with this situation. The systematic procedure for finding a proper 
set of canonical variables for the system has been given by Gardner || whereas Litteljohn 
|H has shown how to construct these variables as power series in an appropriate adiabatic 
parameter 0. The rotational motion of the particle may then be separated from that of 
the guiding center and the guiding center motion result to be described as the phase space 

1 Sij = £ y denotes the completely antisymmetric tensor in two dimensions: en = £22 = 0, £12 = 
—£21 = 1. Everywhere in this paper the sum over repeated indexes is understood. 

2 In his work Lilletjohn made use of a set on noncanonical variables. This may however be reconduced 
to a canonical one by means of an appropriate transformation. 
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motion of the remaining couple of canonical variables. 

We now come to quantum mechanics. The first thing which hits us in the eye is that 
Plank's constants h introduces a length-scale into the problem: the magnetic length 



Quantum mechanics provides therefore an universal means to evaluate the order of mag- 
nitude of the variation of a magnetic field: the comparison between the length-scale over 
which the field varies and the magnetic length. For a rather weak 1 gauss magnetic field, 
e the electron charge, Is is already a very small length, lg ~ 10~ 4 cm. The magnetic 
length Ig appears therefore as a natural adiabatic parameter of the theory, that, on the 
contrary to the classical case, have not to be introduced by hand. 

In this paper we construct a suitable set of operators allowing to rewrite the Hamilto- 
nian describing the motion of a particle in a two-dimensional magnetic field as a formal 
power series in the magnetic length Ib- The expansion parameter having the dimension 
of a length, the perturbative series makes sense only when the n—th derivative of the field 
times Ib" 1 is a very small number in every point of the plane. 

In section 2 we give an outlook to the anomalous Hamiltonian structure of a particle 
interacting with a uniform magnetic field. The construction of a set of operators adapting 
to this simple problem indicates clearly the presence of a fast and a slow degree of freedom 
of the system. It also suggests the way to be followed in the general case. In section 3 we 
construct a new set of variables {IT, X 1 ; i = 1,2} adapting to the inhomogeneous problem. 
In the semiclassical picture the ITs describe the rapid rotation of the particle, whereas 
the X*s take care of the dynamics of the guiding center. The construction is carried on 
order by order in the perturbative parameter Ib up to the second order. In section 4 the 
cinematical momenta 7Tj and the coordinates x l , i = 1,2, of the particle are rewritten in 
terms of the new variables. This allows to write down the Hamiltonian of the system as a 
power series in the magnetic length Ib- An appropriate unitary transformation is finally 
used in section 5 in order to make the perturbative expansion to depend on the variable 
111 and II2 only by means of their combination J = |(IIi 2 + n 2 2 ). The fast rotational 
motion of the particle may so be separated from the slow drift of the guiding center. In 
the classical limit our result reproduces correctly Littlejohn's Hamiltonian ||. Section 6 
contains our conclusions. 




(1) 
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2 An outlook to the canonical structure 



Magnetic interactions appear as modifications of the canonical structure of dynamical 
systems. This feature plays a central role in the Hamiltonian, and hence quantum, de- 
scription of the adiabatic motion of charged particles in strong magnetic fields. We begin, 
therefore, by briefly reviewing this aspect of the problem in the case of a homogeneous 
magnetic field, introducing at the same time the conventions and the notations we use in 
the rest of the paper. 

Let us consider a spin-1/2 particle of mass m, charge e and gyromagnetic factor g 
constrained to move in the x 1 — x 2 plane in presence of the homogeneous magnetic field 
B(af) = Bz. Denoting by cii(x), i = 1,2, an arbitrary choice of vector potential for the 
dimensionless field H(x)/B, diaj — djtii = e^, the dynamics of the particle is described 
by the Pauli Hamiltonian 




« = ^L|-' ? '- "<<-' : > 1 -.'/—"> «■ (2) 

It is customary to parametrize the system by means of the canonical variables pi = 
—ihd/dx 1 and x % , i = 1,2. Nevertheless, Hamiltonian (0) takes an extremely simpler 
form if we replace the canonical momenta Pi by the gauge covariant cinematical momenta 
TTi, TTi = Pi — ^-ai(x). Aside from a scale factor 7Ti and 7C2 behaves as conjugate coordinates 
and Hamiltonian ([|) rewrites as that of an harmonic oscillator, allowing the immediate 
solution of the problem. In order to get an adequate set of canonical variables for the 
description of the system, it is therefore more convenient to keep the cinematical momenta 
7Tj and to replace the coordinates x l by the guiding center operators X\ X % = x % + -^e^tTj. 
In the classical limit the X*s represent the coordinates of the center of the classical orbit. 
The set of operators {TCi,X l ;i = 1,2} fulfills the commutation relations 



TTi, TTj 



X\X J 



MB 

c 

0, (3) 

he 
Jb l 



so that, up to some scale factors, 7Tj and X % may be recognized as canonical variables. We 
note that the cinematical momentum 7i"i and the guiding center coordinates X 2 behave 
now as "canonical coordinates" whereas the cinematical momentum TT2 and the guiding 
center coordinate X 1 correspond to the respective "canonical momenta". The presence 
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of a homogeneous magnetic field produces a kind of rotation of the canonical structure, 
mixing up canonical momenta and coordinates in new canonical variables. 

An outlook to the canonical commutation relations (^) allows to single out a second, 
very important, peculiarity of the variables 7Tj and X 1 . Thinking to Heisenberg's equations 
of motion we immediately realize that the temporal variation of the 7TjS is proportional to 
B, tx, = i[H, Ttil/fi ~ B, whereas that of the X*s goes as 1/B, X 1 = i[H, X^/h ~ In 
the limit of very strong magnetic fields, the canonical structure of the system is distorted 
from the magnetic-scale B and the 7T;S and X l s behave naturally as describing respectively 
a fast and a slow degree of freedom of the system. 

In the next section we construct a set of operators generalizing the 7TjS and the X*s 
for an arbitrary inhomogeneous magnetic field B(x) = Bb(x)z (the dimension of the field 
being reabsorbed in the magnetic-scale B) . This allows to set up an adiabatic description 
of the motion of a quantum particle in a two-dimensional magnetic field in a quite simple 
way. 

For the shake of clearness it results to be convenient to introduce dimensionless quan- 
tities by factorizing the energy scale Hub , ojb = — , from the Hamiltonian. Recalling 

that the magnetic length Is = \f^§, we redefine the cinematical momenta as 

Hamiltonian (§) writes then simply as Ji = %ujb{\ H^f — go's)- It is instead necessary to 
keep guiding center operators with the dimension of a length, 

X i = x l + Ibe^ttj, (5) 

in order to preserve the scale dependence of the canonical commutation relations (|3]) 



7Tj, 71 j 

m,xi 

X\X J 



lEij , 



0, (6) 
-il B 2 e i] . 



A general inhomogeneous magnetic field B(af) depends on coordinates with the dimension 
of a length, after all. 



3 Hamiltonian (Q), of course, does not depend on the X l s. We nevertheless act as it does in order to 
explore the behaviour of the canonical variables in the limit of strong magnetic field. 
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3 Looking for a suitable set of variables 



We now consider a spin-1/2 particle of mass m, charge e and gyromagnetic factor g con- 
strained to move on the x 1 — x 2 plane under the action of the inhomogeneous magnetic 
field B(x) = Bb(x)z. b(x) is supposed to be an arbitrary positive nevervanishing smooth 
function of x. Introducing an arbitrary choice of vector potential a>i(x) for the dimen- 



sionless field b(x)z, diCij 
Hamiltonian 



djdi 



Eijb(x), the system is again described by the Pauli 



(7) 



=1,2 



where the cinematical momenta 7Tj = —ilsdi — aj(x)//s have already been introduced in 
place of the canonical momenta p^. In spite of the very simple dependence of Hamiltonian 
(0) on the 7TjS, the cinematical momenta are no longer conjugate variables and can not be 
directly used to give a simple solution to the problem. The set of operators , x % \ 2 = 1,2} 
fulfills in fact the commutation relations 



7li,X 



x l , X 3 



ieijb(x), 

-il B S J i, 
0. 



(8) 



Nevertheless, we can try to construct a set of variables generalizing the cinematical 
momenta (§) and the guiding center operators @ to the case of an inhomogeneous mag- 
netic field as power series in the magnetic length Is with coefficients depending on the 
7TjS and the x l s. We proceed order by order in the parameter l B by constructing, say at 

2 = 1,2} fulfilling adequate conditions. We 



the n-th order, a set of operators {U.\ ,Xl^ 
require that 

- in the limit of a constant magnetic field, b(x) — > 1, the H^s and the Xl-^s should 

reduce to the cinematical momenta and to the guiding center operators (||D 
respectively, 

— the II^s should be conjugate variables up to terms of order Z# n , 

- the -X"L)S should commutate with the IT^s up to terms of order Ib" '■ 

These commutation relations generalize obviously the @. We nevertheless do not insist 
the XL-jS to be conjugate variables 0]. 
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Zero order variables 



In order to fulfil the conditions above up to terms of order lg we simply rescale the 
cinematical momenta 7Tj by a factor b^ 1 ^ 2 , 



n 



(0) 



(9) 



where the function ft -1 / 2 should be evaluated in X(o) = x. The anticommutator { , } is 
introduced in order to make the Xl\ s hermitian. A brief computation gives the commu- 
tation relation fulfilled by the set of operators {Ilf 1 ,XLy, i = 1,2}, 



n (o) n (o) 

n (0) yj ' 

yi y3 
A (0))^-(0) 



%Eij + lEi 3 
-UnSib- 1 ' 2 , 



(dkb) (o) 
63/2 » u i 



(10) 



0. 



where all the functions have to be evaluated in X^. As required the Ilf^s result to be 

conjugate variables up to terms of order Ib and the commutators between the IQ s and 
the Xfas are again of order Ib- 



First order variables 

We proceed by constructing a couple of operators Xhs, i = 1, 2, commuting with the EQ s 
up to terms of order Ib 2 - The goal is achieved by performing a transformation generalizing 
(||) to the case of an inhomogeneous magnetic field. Xfo is defined as 



4 ) = a r* + ^{6- 1 / 2 ,nW}, 



'11} 



the function b l ^ 2 being evaluated in X(o). We have so the following set of commutation 
relations 



n (o) n (o) 



; 2 



—IE 



'J 





/ (d k b) 


\ b 3 / 2 ' 


km In J 


(d k d t b) 


1 


b 2 



nL 0) + 



b 3 



{n£Mi(°)} \ + o(i B % 



s 



tt(0) yj 



A (i)' A (i) 



il B 2 e jk 



1 ((dtb) 



n 



(0) 



v 2 { b 2 

il B 2 e ij b~ l + o(l B 3 ). 



1 f(gfcg n (o)' 
4 1 6 2 ' 4 



(12) 



All the functions have to be evaluated in X(iy We observe that the XVs are no longer 
commuting variables so that the evaluation of the magnetic field function b and of its 
derivatives in involves ordering ambiguities. The commutator of the Xl^ being of 
order l B 2 , these ambiguities concerns the terms o(l B 3 ). 

The operators 11$ s are constructed by adding a correction of order l B to the lif^s. 
This correction is looked for as a homogeneous second order polynomial in the EQ s. A 
brief computation indicates that the right choice is 



n 



(i) 



nf - -e 
24 



B -km 



(d k b) 

b 3/2 ■ 



{n?\ii£>} 



(13) 



where the functions are evaluated in Xm. The commutation relations fulfilled by the set 



of operators {11$ , -Xm! i = 1, 2} can be immediately obtained by means of ([12]) 



n (1) n (1) 



n (1) yi 



X \l)i X {l) 



(1) ' 

The IT) s behave as conjugate variables up to terms of order l B and the commutators 
between the 11$ s and the Xfas vanish again up to second order terms in l B . 



^B km An 

\l{d k d l b) 5 (0*6) (0,6) r m n , 

i^~ £ £ o— xa q B ' \ U ™, U n j > + °( l B 



il B z e ] 



2 b 2 9 b 3 
3-7* f 1 tt(i)1 1 / ^ rr(i) 



V 2\ b 2 

il B 2 e ij b- 1 + o{l B *). 



n 

4 } b 2 ' 1 



(14) 



Second order variables 

The second and further order variables have to be constructed along the same line of the 
previous one. Looking for the /^-correction to be added to the XLs as a homogeneous 

second order polynomial in the 11$ s we easily find 



x h) 



i 2 

l ^_ £ ik £j i 



b 2 



(15) 
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the functions being evaluated in Xm. The relations ( 1^) allow to compute the commuta- 
tions relations 



n w n (1) " 

LL i '^(2) 

y i yi 

A (2)' A (2) 



tEij ij ^ & ^ 



+ o(l B 3 ), 



\{d k d x h) 5 (0*6)® 6) 



2 6 2 
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fc 3 



,{n« n«} 



(16) 



where every function is now evaluated in X{n\. Finally we proceed by constructing II 
by adding an adequate homogeneous third order polynomial in the IT-^s to , 



(2) 



n 



(2) 



n (D + ^ £ km £ ln 



l {d k dib) 5 (d k b)(dib) (1) ( i) 
2 b 2 9 b 3 ' m i n 



(17) 



the function b and its derivatives being again evaluated in Xp)- The set of operators 



{nfU, 



(2)' 



1, 2} fulfills the desired commutation relations 

ieij + o(/ B 3 ), 



n (2) n (2) 



n ( 2 ) yi 
yi yi 

A (2)>^-(2) 



+ o(l B 3 ), 

-ilB^b' 1 +0(l B 3 ). 



(2) ■ 

The II- s are conjugate up to terms of order l B and commutate with the Xfys again 
up to the third order in l B . In the limit of a constant magnetic filed, of course, all the 
variables that we have introduced have the correct behaviour. 



A noncanonical set of variables 

In principle we may think to repeat this procedure an arbitrary number of times and to 
construct, as power series in l B with coefficients polynomial in the 7TjS and depending on 
the x l s through the function b and its derivatives, a set noncanonical operators {IT, X 1 ; i = 
1,2} fulfilling the commutation relations 



IIi> YLj XEij, 



IT,X J 



X\X J =-il B 2 e ij b~ 



(19) 
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b^ 1 being now evaluated in X. Anyway, in order to discuss the problem to the second 
order in the adiabatic parameter Is, we only need to know the first three terms of these 
series, terms which we have already evaluated. Up to terms of order /# 3 we may therefore 
confuse the EE; s and XL-.S with the IljS and X l s respectively. 

Let us observe now, that aside from the function 6 _1 appearing in the right-hand side of 
the third identity, the commutation relations (|19D correspond to the Although the X l s 
are not conjugate variables they commute with the Us and their commutator is a function 
of the A*s alone. The Hilbert space of the system, separates therefore under the action 
of the two couple of operators in the direct sum of two subspaces describing each one a 
degree of freedom of the system. The scale dependence of the commutation relations ([19|) 
indicates again that in the limit of very strong magnetic fields (small Ib) the IljS and the 
X l s describe respectively a fast and a slow degree of freedom of the system. The adiabatic 
cinematical momenta Hi and the adiabatic guiding center operators X 1 introduced in this 
section, appear therefore as a suitable choice of variables for our problem. 

4 The adiabatic expansion 

The next steep is to rewrite Hamiltonian (^]) in terms of the new set of variables. To 
this purpose we have first to invert the power series expressing the IljS and the X*s in 
terms of the 7TjS and the x l s and then to replace these expressions in (0). As a result the 
Hamiltonian will appear as a power series in the magnetic length Ib 

n = n (0) + i B n {1) + i B 2 n {2) + ... . (20) 

In the limit of slowly-varying magnetic field, that is when the magnetic length may be 
considered small with respect to the length-scale over which the magnetic fields varies, 
this equation may be interpreted as a perturbative expansion of the hamiltonian and used 
to get approximate expressions of the spectrum and the wavefunctions of the system. We 
call it the adiabatic expansion of the Hamiltonian. 

The task of obtaining the operators 7Tj and x l in terms of the new variables Hj and X % 
is not a hard one. Recalling that X 1 = XL-, + o(Ib 3 ) and using equations (|9|), (|TTD, (|T3| ) 
and (0) we can immediately obtain the first three terms of the power series expressing 
the x l s in terms of the 7TjS and the X*s, x l = x 1 (tt, X). This allows to rewrite the n^s as 
functions of the the 7TjS and the X l s, so that by using equations QT^D, ( |TT| ) and the relation 

f 2) 

= II, • + o(Ib ), is possible to rewrite the n^s as power series in Ib with coefficients 
depending on the 7TjS and the X z s, Hi = EL, (7?, X). Inverting these series order by order 
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we can finally calculate the 7TjS as functions of the IljS and the X*s. By substituting these 
expressions in x l = x l (7T,X), we also get the x l s as functions of the IL„s and the X l s. A 
few computations conduce to the result 

7* = 6 1 / 2 n l -|M^ {n ,,n m } + 

x l = X^-he^b-^-^eV^iU^ + oih"), (22) 

12 o z 

the function b and its derivatives being evaluated in X. As a check we can easily reobtain 
the commutation relations (]8|) by means of the expressions (J2"T|), (^) and (|Hf). 

The first three terms of the adiabatic expansion of the Hamiltonian can now be im- 
mediately obtained by substituting the expressions (^l|) and (|22|) in (^). In order to 



simplify the notation we introduce the operators J = I^Tljllj, = -S^TlilliJIj and 
Jki — ^S^Hi {Ilfc, Ih} IIj. J is just the Hamiltonian of an harmonic oscillator in the con- 
jugate variables IT!, n 2 and, as we will see in the next section, plays a central role in our 
discussion. A brief computation yields 

H i0) /huj B = b (J + ga 3 ) (23) 
H^/huB = + (24) 



1 (didjb) 5 fab) (djb)' 



16 b 144 b 2 



S ij (25) 



These equations represent a good starting point for a perturbative analysis of the system. 
Nevertheless there is still some more work that may be done in order to make this task 
simpler. In the next section we will follow a strategy which very much recalls the one 
used in the perturbation theory of Hamiltonian mechanics, in order to find a unitary 
transformation making Tt to depend on Hi and U2 only by means of their combination J 
and its powers. J results then into an adiabatic invariant which may be identified with 
magnetic moment of gyration of the particle. 
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5 Effective guiding center dynamics 



A well known strategy in treating perturbative problems in classical mechanics is that 
of subjecting the the system to a series of near-identity canonical transformations in 
order to make the various orders of the perturbative expansion of the Hamiltonian to 
be independent on certain variables. In this way is possible to individuate the constants 
of motion of the system up to an arbitrary order in the perturbative parameter. The 
approximate solution of the problem results then notably simplified. In this section we 
will adopt the equivalent technique in quantum mechanics. Hamiltonian ( POD will be 
subjected to a near-identity unitary transformation such that the following criteria is 
fulfilled: the new Hamiltonian should depend on the conjugate variables Hi and H 2 only 
by means of the operator J and its powers. In this way J results into a constant of motion 
of the system and the dynamics of the X l s separates from that of the IljS. 

Let us therefore consider a near-identity unitary operator U, represented in the form 
U = e lC . The hermitian operator C will differ from the identity operator 1 by terms of 
order grater than lg so that we may represent it as the power series 

l + l B £W+l B *£W + ... . 



c 

Subjecting Hamiltonian (p0|) to this transformation produces the result 



(26) 



+i B 2 n^ + i\c {2 \n^ 



(27) 



The zero order term of the Hamiltonian remains unchanged whereas the others are cor- 
rected by additive terms depending on the commutators of the & n >$ with and the 
subsequent terms of the adiabatic expansion. We start therefore by looking for an op- 
erator £W such that its commutator with H' ', [& 1 \Tv-°>], might correct the first order 
term of the adiabatic expansion in the desired way. A brief computation indicates that it 
is possible to annihilate the first order term of the transformed Hamiltonian by choosing 



£(1) = _ 



{9S „ 7 , 



&V2 



<P* (2J k /3 + ga 3 U k ) 



The next steep is to construct the operator making the second order term of the 
adiabatic expansion, H {2) + i[C^\H^] + i[C^\ H {1) ]/2, to depend on n 1; n 2 , only by 
means of their combination J and its powers. The correct choice results to be 



£(2) 



'(didjb) lidMdjb)' 



16 



b 2 



b 3 



e ik V l {J kl + 2ga 3 {n fc ,nj) 



(29) 
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Having obtained further terms of the adiabatic expansion it would be possible to 
keep on this procedure to an arbitrary order. At least in principle therefore is possible 
to make J into a constant of motion to an arbitrary order in the adiabatic parameter 
Ib- Our analysis stops to the second order. Subjecting Hamiltonian (|20|), with TiS°\ 
and Ti.^ given respectively by ( p3f) , ([24] ) and (p5j ), to the unitary transformation already 
described we get the new adiabatic expansion 



H'/Hub = b (J + ga 3 ) + 



/ 2 



A6 |V6| 2 



+■ 



7 2 r 
JL 
16 



(j 2 + 2^3 j) + 



6 2 



+ o(/ B 3 ) 



(30) 



where the functions b, Ab = 5 % ^{didjb) and | V6| 2 = 5 % ^{dib){djb) are again evaluated in X. 
Up to terms of order Ib 3 the operator J results into a constant of motion of the system 
and the harmonic oscillator degree of freedom described by Hi and n 2 separates from the 
the guiding center motion described by X 1 and X 2 . In the semiclassical picture the fast 
rotation of the particle around the guiding center corresponds to the phase space motion 
in the plane 111 — H2 whereas the slow drift of the guiding center in the plane corresponds 
the phase space motion in the plane X 1 — X 2 . The adiabatic invariant J may therefore 
be identified with the magnetic moment of gyration of the particle, that is the magnetic 
moment of the current loop described by the particle in one cyclotron gyration. In the 
classical limit the Hamiltonian operator (^) reproduces correctly Littlejohn's guiding 
center Hamiltonian H. 

Once the system has been frozen in one of its gyrating eigenstates, that is, in an 
eigenstate of the fast degree of freedom, Hamiltonian ( ftCD describes the corresponding 
effective guiding center dynamics. It is interesting to observe that the zero order effective 
guiding center Hamiltonian corresponds to the magnetic field function b(x) in which the 
coordinates x l and x 2 have been substituted with the noncommuting operators X 1 and 
X 2 . Introducing a pair of Euler potentials x x (x) and x 2 (x) for the magnetic field b(x)z, 
Vx 1 AVx 2 = b P|, it is also possible to make the X l s into a couple of conjugate operators. 



Defining X 1 = x\X) and X 2 = x 2 (X) we get in fact [X\X 2 ] = -il B 2 , g |. In order 
to evaluate the zero order spectrum of the system it is therefore sufficient to construct 
the Hamiltonian operator of the system by substituting the commuting variables x 1 and 
x 2 with the couple of conjugate operators X 1 and X 2 in the function ^(^(x^x 2 )). This is 
similar to a quantization procedure, which for some aspects, reminds the one explored by 
J. R. Klauder 1T2L IXBl. 
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Further corrections to the spectrum may be evaluated by means of standard pertur- 
bation theory and the subsequent terms of the adiabatic expansion (j3C|). 

To conclude the section we specialize to the case of electrons, g = —1/2, in order to 
check our result to be consistent with a theorem of Y. Aharonov and A. Casher stating, 
among other things, that the ground state of an electron moving in an arbitrary two- 
dimensional magnetic field has always zero energy ||. By setting the gyromagnetic factor 
to —1/2 our Hamiltonian reduces to the quite compact form 

H'/hujB = b (j- i ff3 ) + ^ 

Recalling that the eigenvalues of J are E n = (n + |), n = 0, 1,2..., and that of cr 3 are 
s = ±1 it is immediate to observe that the ground sate energy of Hamiltonian (|3lD , 
(n = 0, s — 1), is zero up to terms of order l B 3 . 



Ab JV&I 2 



J - er 3 J + 



+ o(l B 3 ) 



(31) 



6 Discussion and conclusion 

The construction of a set of noncanonical operators sharing the main properties of the 
cinematical momenta and the guiding center operators of a particle moving in a homoge- 
neous magnetic field, allows to separate the unperturbed dynamics of a quantum particle 
moving in a two-dimensional weakly-inhomogeneous magnetic field, from the perturbation 
making the guiding center dynamics nontrivial. The effective guiding center Hamiltonian 
writes furthermore as a power series in the magnetic length l B . The main features of the 
construction are that it is perturbative in nature, that is, it is performed order by order 
in the adiabatic parameter l B , and involves only simple algebraic manipulations. These 
properties allows to achieve the goal in a very simple and economical way. 

It is worth-while to spend some words on the efficiency of the adiabatic expansion 



(30). As we said in the introduction the presence of a magnetic field in a quantum context 
introduces the length-scale l B . Intuitively, we may therefore consider a magnetic field to be 
slowly- varying (weakly inhomogeneous) , when its variation ratio over the magnetic length- 
scale l B is a small number. The adiabatic expansion obtained in this paper completely 
confirms this picture. Furthermore , still for a rather week magnetic field of 1 gauss, the 
magnetic length is very small, l B ~ 10 cm. To have one's eye on some number, we recall 
that the geomagnetic field near to the surface of the earth is of the order of the gauss as 
well as the magnetic field produced by a wire in which flows a current of a few ampere. 
In the physics of the quantum Hall effect much stronger fields are used, typically of the 
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order of 10 4 — 10 5 gauss. All these fields varies significantly over lengths which go from 
the millimeter to some meters. That is to say, the adiabatic approximation is a very good 
one. Among the other things, therefore, our work gives a further explanation of the why 
so sharp Landau levels are observed in the quantum Hall effect. 

Once the system has been frozen in one of its gyrating eigenstates, the dynamics of 
the remaining degree of freedom is described by means of a Hamiltonian operator which 
have to be constructed by quantizing the classical Hamiltonian 6(af(x 1 ,x 2 )) (replacing h 
with Ib )■ The Ib -terms of the adiabatic expansion ( p0|) will in general contribute a small 
correction to the spectrum of the system. By an appropriate choice of the magnetic field, 
therefore, an arbitrary one-degree-of-freedom Hamiltonian may be reproduced. It would 
be interesting to explore the possibility of concretely realize devices of this kind. 

We conclude by observing that the method we developed in this paper works in the full 
three-dimensional well as in the two-dimensional one. The study of the adiabatic 

motion of a charged, spinning, quantum particle in a three-dimensional magnetic field 
will be reported upon in a forthcoming publication. 
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